Soliton-induced optical absorption of halogen-bridged mixed-valence 
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Employing the one-dimensional single-band extended Peierls-Hubbard model, we investigate 
optical conductivity for solitonic excitations in halogen-bridged binuclear metal (MMX) complexes. 
Photoinduced soliton absorption spectra for MMX chains possibly split into two bands, forming 
a striking contrast to those for conventional mononuclear metal (MX) analogs, due to the broken 
electron-hole symmetry combined with relevant Coulomb and/or electron-phonon interactions. 



PACS numbers: 71.45.Lr, 42.65. Tg, 78.20.Ci, 78.20.Bh 
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I. INTRODUCTION 



Halogen (X)-bridged mixed-valence metal (M) com- 
plexes [1,2] such as [Pt(en)2X](C104)2 (X = CI, Br, I; 
en = ethylenediamine = C2H8N2), which are re- 
ferred to as MX chains, have been playing a prominent 
role in understanding the electronic properties of one- 
dimensional Peierls-Hubbard systems. The competing 
electron-electron and electron-phonon interactions yield 
various ground states [3,4], which can be tuned by chemi- 
cal substitution [5] and pressure [6]. Solitonic excitations 
inherent in charge-density-wave (CDW) ground states 
stimulate further interest in AIX complexes. In this con- 
text, we may be reminded of polyacetylene, the trans 
isomer of which exhibits topological solitons [7,8]. Sev- 
eral authors [9-11] had an idea of similar defect states 
existing in MX chains. Photogenerated solitons [12,13] 
were indeed observed in the PtX compounds. 

In recent years, binuclear metal analogs of MX com- 
plexes, which are referred to as MMX chains, have at- 
tracted further interest exhibiting a wider variety of 
ground states [14,15], successive thermal phase transi- 
tions [16,17], photo- and pressure-induced phase transi- 
tions [18-20], and incomparably larger room-temperature 
conductivity [21]. In such circumstances, soliton so- 
lutions of an MMX Hamiltonian of the Su-Schrieffer- 
Heeger type [7] have recently been investigated both an- 
alytically and numerically [22]. The direct M~M over- 
lap contributes to the reduction of the effective on-site 
Coulomb repulsion and therefore electrons can be more 
itinerant in MMX chains. Hence we take more and more 
interest in solitons as charge or spin carriers. 

The ground-state properties of MMX complexes were 
well revealed by means of X-ray diffraction [23] , nuclear 
magnetic resonance [24] , and the Raman and Mossbauer 
spectroscopy [16], while very little [25,26] is known about 
their excitation mechanism. Photoinduced absorption 
spectra, which served as prominent probes for nonlin- 
ear excitations in MX complexes [12,13], have, to our 
knowledge, not yet measured on MMX complexes prob- 
ably due to the lack of a guiding theory. Thus moti- 
vated, we study optical conductivity for MMX solitons 



with particular emphasis on a contrast between topical 
MMX and conventional MX complexes. Photoexcited 
MMX chains may yield distinct spectral shapes due to 
the definite breakdown of the electron-hole symmetry. 



II. MODEL HAMILTONIANS AND THEIR 
GROUND-STATE PROPERTIES 



We describe MX and MMX chains by the one- 
dimensional i- and |-filled single-band Peierls-Hubbard 
Hamiltonians 
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Here, n^-n^s = ^'^^_i:n,s'^ti-n,s with ajj .„ being the creation 
operator of an electron with spin s — ±^ (up and down) 
for the M d^^ orbital labeled as /i = 1, 2 in the nth MX 
or MMX unit, Imm and Imxm describe the intra- and 
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(a) CDW in MX (b) CDW in MMX 

FIG. 1. Schematic representation of CDW ground states 
observed in MX (a) and MAIX (b) complexes. 
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complexes, we enforce the constraint 
every M — X — M bond. Then the force equilibrium con- 
ditions 



dl 
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interunit electron hoppings, respectively, a and /3 are the 
site-ofF-diagonal and site-diagonal electron-lattice cou- 
pling constants, respectively, In^ — Vn — Un-i and 
li^"^ = Un — Vn with Un and Vn being, respectively, the 
chain-direction displacements of the halogen and metal in 
the nth unit from their equilibrium position, and Kmx 
is the metal-halogen spring constant. We assume, based 
on the thus-far reported experimental observations, that 
every M2 moiety is not deformed. The notation is further 
explained in Fig. 1. We set tMXM and Kmx both equal 
to unity in the following. 

The existent MMX complexes consist of two families: 
i?4[Pt2(pop)4X]-nH20 (X = Cl,Br,I; i? = Li, K, Cs, • • •; 
pop = diphosphonatc = P2 05H|^) and M2(dta)4l 
(M = Pt,Ni; dta = dithioacetate = CH3CS2~), which 
are referred to as pop and dta complexes. The pop com- 
plexes possess the mixed-valence ground state illustrated 
in Fig. 1(b) [24,27], while the dta complexes exhibit dis- 
tinct types of ground states [14-16,28] due to the pre- 
dominant site-off-diagonal electron-phonon coupling [29] 
or metal-on-site Coulomb repulsion [30]. Conventional 
MX complexes structurally resemble MMX complexes 
of the pop type and have the mixed-valence ground state 
illustrated in Fig. 1 (a) [2] , which is symmetrically equiv- 
alent to that in Fig. 1(b) [15,31]. Hence we study the 
pop complexes in our first attempt to compare the photo- 
products in MMX and MX chains. While the dta com- 
plexes behave as d-p-hybridized two-band materials in 
general, the pop complexes are well describable within 

-single-band Hamiltonians [15]. Since the site-off- 
diagonal electron-phonon coupling is of little significance 
with the CDW backgrounds of the Fig. 1 type, we set a 
equal to zero. The rest of parameters are, unless other- 
wise noted, taken as Um — 1-2, Vmxm = 0.3, and [3 = 0.7 
for MX chains [1,32,33], while Um = 1-0, Vmm = 0.5, 
Vmxm = 0.3, and /3 = 1.4 for MMX chains [19,26]. 

III. CALCULATIONAL PROCEDURE 

We treat the Hamiltonians (1) within the Hartree-Fock 
approximation. The lattice distortion is adiabatically de- 
termined and can therefore be expressed in terms of the 
electronic density matrices {n^;n,s)BFi where (■ • •)hf de- 
notes the thermal average over the Hartree-Fock eigen- 
states. Since no spontaneous deformation of the metal 
sublattice has been observed in any MX and MMX pop 



"^KmxIu ^ = /3y^((n^i:K,s)HF - {nu:n-l.s)wF), 
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^KMxln'^ = /3^((n-i/:n,s)HF — (n^:„+l,s)HF), 



where 11 = v = 1 for MX chains, while /i = 1, v = 2 for 
MMX chains. 

The real part of the optical conductivity is given by 
^(^) = ]^ E - fie')W\J\^)\'S{e' - e - hu;), 

e,e' 

(8) 

where /(e) = (e'/'^^^^ + l)-i and {e'\J\e) is the matrix 
element of the current density operator J between the 
eigenstates of energy e and e'. J' is defined as 

Jmx = Aixli-^ Jmmx = Aim + Jmxm^ (9) 
for MX and MMX chains, respectively, with 

3 MM = -^CMMtMM y^ial:n.sO-t^-n,s - al:n,sO-i^:n,s) , (10) 
n,s 

JmXM = J^CMXM ^ [tMXM - "(4+1 + 4^^)] 
n.s 

where cmm and cmxm are the average M — M and 
M—X—M distances, respectively, and are set for cmxm — 
2cMM [34,35]. 

In order to elucidate the intrinsic excitation mecha- 
nism, solitons are calculated at a sufficiently low temper- 
ature without any assumption on their spatial configu- 
rations. Charged solitons (S'^; a — ±) are obtained by 
setting the numbers of up- and down-spin electrons, A^-|_ 
and both equal to {p — 1/2)N — a/2, while spin-s 
neutral solitons [S°*; s = ±i] with N± = {p-l/2)N±s, 
where N, taken to be 401 in our calculation, is the num- 
ber of unit cells and p is the number of metal sites in the 
unit cell. 
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FIG. 2. Spatial configurations of tlie optimum solitons in 
MX and MMX chains. 
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FIG. 3. Tfie formation energies of the stably located soli- 
tons as functions of the Coulomb interaction, where /3 = 0.7 
and Um varies keeping the relation Um = Vmxm /0.25 for 
MX chains, while /3 = 1.4 and Um varies keeping the relation 
Um = Vmm/O.S = Vmxm /0.3 for MMX chains. 



IV. RESULTS 

We show the spatial configurations of the optimum 
solitons in Fig. 2 and plot their formation energies AEs 
in Fig. 3. In MX chains, charged solitons are stable 
laying their centers on halogen sites, while neutral soli- 
tons on metal sites. In MMX chains, all the optimum 
solitons are halogen-centered. As the Peierls gap Sgap in- 
creases, solitons generally possess increasing energies and 
decreasing extents, and end up with immobile defects. 
Um and Vmm reduce i?gap and thus enhance the mobil- 
ity of solitons, while Vmxm has an opposite effect. With- 
out any Coulomb interaction, the soliton formation en- 
ergies are all degenerate and scaled by the Peierls gap as 
AEs = Egs^p/ir in the weak-coupling region i?gap ^ tMXM 
[36], whereas their degeneracy is lifted and neutral soli- 
tons have higher energies than charged solitons with fur- 
ther increasing gap. Since Sgap decreases with increasing 
Um, AEs is a decreasing function of Um- Neutral soli- 



tons are more sensitive to the Coulomb interaction and 
thus their formation energy turns smaller than that of 
charged solitons with increasing Um- Considering that 
Sgap is a linear function of Um , we learn that neutral soli- 
tons well keep the scaling relation AEs oc i?gap against 
the Coulomb interaction. 

Figure 2 suggests that an MX soliton of charge a and 
spin s, S'^'*, described in terms of electrons is equivalent to 
its counterpart S"*^"** described in terms of holes. Such a 
symmetry is more directly observed through the energy 
structures shown in Fig. 4. Solitons generally exhibit 
an additional level within the gap. There appear further 
soliton-related levels in the strong-coupling region [37]. 
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FIG. 4. Density of states for the optimum soliton solutions. 
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FIG. 5. (Color online) Energy shifts of the localized soli- 
ton levels from the gap center scaled by the Peierls gap as 
functions of the Coulomb interaction. The parametrization is 
the same as that in Fig. 3. The level structures of up- and 
down-spin electrons are plotted by solid (broken) and broken 
(solid) lines, respectively, for S"^ (S"^), while they are degen- 
erate for charged solitons. 
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FIG. 6. Absorption spectra with a soliton excited for MX 
and MMX chains. The low-frequency structures are scaled 
up in insets. 



The intragap soliton levels are analyzed in more detail 
in Fig. 5. When the coupling strength increases without 
any Coulomb interaction, the mid-gap level due to a neu- 
tral soliton keeps still, while the charged-soliton levels de- 
viate from the gap center, symmetry. Once the Coulomb 
interaction is switched on, the neutral-soliton level also 
begins to move away from the gap center breaking the 
spin up-down Soliton-related electron and hole levels are 
necessarily symmetric with respect to the center of the 
gap in MX chains, while such a symmetry is generally 
broken in MMX chains. The electron- hole symmetry is 
kept and broken in the MX Hamiltonian (la) and the 
MMX Hamiltonian (lb), respectively. 

Now we show in Figure 6 photoinduced absorption 
spectra which arc expected of MX and MMX chains. 
Since photoinduced solitons are necessarily in pairs, we 
superpose the (S°^) spectrum on the S+ (S°^) spec- 
trum. The spectra of S+ and are degenerate for 
MX chains but distinguishable for MMX chains. The 
asymmetry between the absorption bands can remain 
even in the absence of the Coulomb correlation provided 
the electron- phonon interaction is sufficiently strong. Al- 
though the spectra of 8°^ and S^^ are degenerate for both 
MX and MMX chains, those for MMX chains are inter- 
esting in themselves. The S'^ spectrum is single-peaked 
for MX chains but splits into two definite peaks for MMX 
chains. The doublet structure of the S° spectrum is due 
to the combination of the broken electron-hole symme- 
try and relevant Coulomb correlation. In principle MX 
chains also lose the electron-hole symmetry with their X 
Pz electrons activated [38], which indeed applies to the 
NiX compounds [39]. However, it is not the case with the 
VtX and PdX compounds, where the energy level of the 
X Pz orbitals is much lower than that of the M d^^ or- 
bitals and thus the p-orbital contribution may effectively 
be incorporated into the intermetal supertransfer energy 
of the single-band Hamiltonian (la). In fact a two-band 



calculation of the PtCl compound [40] still results in a 
single-peaked S° absorption band near the gap center. 



V. SUMMARY AND DISCUSSION 



Since polarons are generated only from electron-hole 

pairs with large excess energies, lower-energy-excited 
charge transfer excitons relax into soliton pairs or de- 
cay by luminescence [2]. A neutral-soliton pair has 
a lower formation energy [1,22,33,41] but a shorter 
life time [2] under realistic parametrizations. Pho- 
toexcited [Pt(en)2X](C104)2 {X = CI, Br) [42,43] and 
[Pt(en)2l](C104)2 [44] indeed yield intragap absorptions 
attributable to neutral and charged solitons, respectively. 
The charged-soliton spectrum has a singic-pcakcd Gaus- 
sian shape, which is well consistent with Fig. 6, while the 
neutral-soliton spectrum consists of a midgap main band 
and an accompanying shoulder structure, which demands 
further driving forces such as lattice fluctuations [32]. 

Solitons may be photogenerated in MMX complexes 
as well, yielding a contrastive absorption spectrum of 
doublet structure. The separation between the two S° 
bands is in proportion to the on-site Coulomb repul- 
sion and may therefore be reduced with pressure applied, 
whereas the asymmetry between the two S^ bands is 
sensitive to the electron-lattice coupling and may thus 
be enhanced by the halogen replacement CI — > Br — > I. 
Photoinduced absorption measurements on MMX com- 
plexes may not only reveal the novel doublet structures of 
soliton-induced spectra but also give a key to the unset- 
tled problem in MX complexes— the shoulder structure 
of the S'' absorption spectrum. Let us make a close col- 
laboration between experimental and theoretical investi- 
gations of dynamic properties of MMX complexes. 



ACKNOWLEDGMENTS 



The authors are grateful to K. Iwano and H. Okamoto 
for fruitful discussions and helpful comments. This work 
was supported by the Ministry of Education, Culture, 
Sports, Science, and Technology of Japan, the Nissan 
Science Foundation, and the Iketani Science and Tech- 
nology Foundation. 



[1] J. T. Gammel, A. Saxena, I. Batistic, A. R. Bishop, and 
S. R. Phillpot, Phys. Rev. B 45, 6408 (1992). 

[2] H. Okamoto and M. Yamashita, Bull. Chem. Soc. Jpn. 
71, 2023 (1998). 

[3] K. Nasu, J. Phys. Soc. Jpn. 52, 3865 (1983). 



4 



I. Batistic, J. T. Gammol, and A. R. Bishop, Phys. Rev. [25 
B 44, 13228 (1991); H. Roder, A. R. Bishop, and J. T. 
Gammel, Phys. Rev. Lett. 70, 3498 (1993). [26 
H. Okamoto, T. Mitani, K. Toriumi, and M. Yamashita, 
Mater. Sci. Eng. B 13, L9 (1992). [27 
J. T. Gammel and G. S. Kanner, Synth. Met. 70, 1191 
(1995). 

W. P. Su, J. R. Schrioffcr, and A. J. Heeger, Phys. Rev. 

Lett. 42, 1698 (1979); Phys. Rev. B 22, 2099 (1980). [28 

H. Takayama, Y. R. Lin-Liu, and K. Maki, Phys. Rev. B 

21, 2388 (1980). [29 

S. Ichinose: Sohd State Commun. 50 (1984) 137. 

Y. Onodera: J. Phys. Soc. Jpn. 56 (1987) 250. [30 

D. Baeriswyl and A. R. Bishop, J. Phys. C 21, 339 (1988). 

N. Kuroda, M. Sakai, Y. Nishina, M. Tanaka, and S. [31 

Kurita, Phys. Rev. Lett. 58, 2122 (1987). 

H. Okamoto, T. Mitani, K. Toriumi and, M. Yamashita, [32 
Phys. Rev. Lett. 69, 2248 (1992). 

M. Kuwabara and K. Yonomitsu, Mol. Cryst. Liq. Cryst. [33 

341, 533 (2000); 343, 47 (2000); J. Phys. Chom. Sohds 

62, 435 (2001). [34 

S. Yamamoto, Phys. Lett. A 258, 183 (1999); J. Phys. 

Soc. Jpn. 69, 13 (2000); Phys. Rev. B 63, 125124 (2001). [35 

H. Kitagawa, N. Onodera, T. Sonoyama, M. Yamamoto, 

T. Fukawa, T. Mitani, M. Scto, and Y. Macda, J. Am. [36 

Chcm. Soc. 121, 10068 (1999). [37 

S. Yamamoto, J. Phys. Soc. Jpn. 70, 1198 (2001). 

B. I. Swanson, M. A. Stroud, S. D. Conradson and M. [38 

H. Zietlow, Sohd State Commun. 65, 1405 (1988); M. A. 

Stroud, H. G. Drickamcr, M. H. Zietlow, H. B. Gray, and [39 

B. I. Swanson, J. Am. Chem. Soc. Ill, 66 (1989). 

S. Yamamoto, Phys. Rev. B 64, 140102(R) (2001). 

H. Matsuzaki, T. Matsuoka, H. Kishida, K. Takizawa, H. 

Miyasaka, K. Sugiura, M. Yamashita, and H. Okamoto, [40 

Phys. Rev. Lett. 90, 046401 (2003). 

H. Kitagawa, N. Onodera, J.-S. Ahn, and T. Mitani, [41 
Synth. Met. 86, 1931 (1997). [42 
S. Yamamoto and M. Ichioka, J. Phys. Soc. Jpn. 71, 189 
(2002); S. Yamamoto, Mol. Cryst. Liq. Cryst. 379, 555 [43 
(2002). 

S. Jin, T. Ito, K. Toriumi, and M. Yamashita, Acta Cryst. 

C 45. 1415 (1989). [44; 
N. Kimura, H. Ohki, R. Ikeda, M. Yamashita, Chem. 
Phys. Lett. 220, 40 (1994). 



Y. Wada, T. Furuta, M. Yamashita, and K Toriumi, 
Synth. Mot. 70, 1195 (1995). 

M. Kuwabara and K. Yonemitsu, J. Mater. Chem. 11, 
2163 (2001). 

L. G. Butler, M. H. Zietlow, C.-M. Che, W. P. Schae- 
fer, S. Sridhar, P. J. Grunthaner, B. I. Swanson, R. J. 
H. Clark, and H. B. Gray, J. Am. Chcm. Soc. 110, 1155 
(1988). 

S. A. Borshch, K. Prassides, V. Robert, and A. O. Solo- 
nenko, J. Chem. Phys. 109, 4562 (1998). 
C. Bellitto, A. Flamini, L. Gastaldi, and L. Scaramuzza, 
Inorg. Chem. 22, 444 (1983). 

C. Bellitto, G. Dessy, and V. Fares, Inorg. Chem. 24, 
2815 (1985). 

S. Yamamoto, Phys. Lett. A 247, 422 (1998); Synth. 
Met. 103, 2683 (1999). 

K. Iwano and K. Nasu, J. Phys. Soc. Jpn. 61, 1380 
(1992). 

Y. Tagawa and N. Suzuki, J. Phys. Soc. Jpn. 64, 2212 

(1995) . 

C.-M. Che, F. H. Herbstein, W. P. Schaefer, R. E. Marsh, 

and H. B. Gray, J. Am. Chem. Soc. 105, 4604 (1983). 

R. J. H. Clark, M. Kurmoo, H. M. Dawes, and M. B. 

Hursthousc, Inorg. Chcm. 25, 409 (1986). 

S. Yamamoto, Phys. Rev. B 66, 165113 (2002). 

Y. Tagawa and N. Suzuki, J. Phys. Soc. Jpn. 59, 4074 

(1990). 

J. T. Gammel, R. J. Donohoe, A. R. Bishop, and B. I. 

Swanson, Phys. Rev. B 42, 10566 (1990). 
H. Okamoto, Y. Shimada, Y. Oka, A. Chainani, T. Taka- 
hashi, H. Kitagawa, T. Mitani, K. Toriumi, K. Inoue, 
T. Manabe, and M. Yamashita, Phys. Rev. B 54, 8438 

(1996) . 

S. W. Wcbcr-Milbrodt, J. T. Gammel, A. R. Bishop, and 

E. Y. Loh, Jr., Phys. Rev. B 45, 6435 (1992). 

K. Iwano, J. Phys. Soc. Jpn. 66, 1088 (1997). 

N. Kuroda, M. Ito, Y. Nishina, and M. Yamashita, J. 

Phys. Soc. Jpn. 62, 2237 (1993). 

H. Okamoto, Y. Kaga, Y. Shimada, Y. Oka, Y. Iwasa, 
T. Mitani, and M. Yamashita, Phys. Rev. Lett. 80, 861 
(1998). 

H. Okamoto, Y. Oka, T. Mitani, and M. Yamashita, 
Phys. Rev. B 55, 6330 (1997). 



5 



